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We demonstrate experimentally that disorder enhanced Andreev current in a tunnel junction between a normal 
metal and a superconductor provides a method to measure electronic temperature, specifically at temperatures 
below 200 mK when aluminum is used. This Andreev thermometer has some advantages over conventional 
quasiparticle thermometers: for instance, it does not conduct heat and its reading does not saturate until at 
lower temperatures. Another merit is that the responsivity is constant over a wide temperature range. 


Thermometers are a cornerstone in experimental 
physics, from the premises of the thermodynamics to the 
latest observations of the universe involving bolometer^. 
Yet in experiments in the millikelvin range and for ob¬ 
jects of nanometer scales, it is hard to measure the actual 
temperatur^. Indeed, in mesoscopic transport measure¬ 
ment, local temperature variations can exist, self-heating 
needs to be avoided and the electronic temperature might 
deviate from the bath temperature. This leads to a need 
of special kinds of thermometers taking such constraints 
into consideration. 

Tunnel junctions have been used in this context for a 
long tim^, both with DC and more recently with fast 
readout!^, but they suffer from various limitations. Op¬ 
erated at finite bias they heat (or cool) the circuit in 
which they are embedded iiP. In addition, their respon¬ 
sivity tends to saturate at low temperature due to par¬ 
ticular tunneling processes, such as environment assisted 
tunneling or Andreev reflection^ . 

We are proposing to use the disorder enhanced An¬ 
dreev current as a temperature probe. When a single 
electron attempts to leave the normal metal (N) trough 
a tunnel barrier (I), it is usually reflected if no states are 
available in the superconductor (S), i.e. if eY < A. Due 
to elastic scattering, this reflected electron may bounce 
around and finally make one or several attempts to tun¬ 
nel. If the coherence time is long enough, all these at¬ 
tempts sum up coherently, and the rare case of being re¬ 
flected as a hole (emitting a Cooper pair in S), becomes 
probable. This leads to a finite current, even for eV < A. 
Measuring this current should give us a direct measure 
of the electronic temperature in the normal metal. The 
response is expected to be a linear function of temper¬ 
ature, which is an advantage, but the main reward is 
its operation point near zero bias. Therefore, such ther¬ 
mometer does not produce much self-heating, so one can 
easily think to use it in measurements where even tiny 
back-action from the thermometer to the system is detri¬ 
mental. To demonstrate the principle of a thermometer 
based on Andreev current, we fabricated a sample de- 
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vice using electron beam lithography and shadow angle 
deposition. The junctions are formed by in-situ oxida¬ 
tion of the Hs = nm thick aluminum leads, before 
evaporating the island in the same chamber. The island 
is based on a titanium gold bilayer, so that the full is¬ 
land behaves as a normal metal. An aluminum/titanium 
buffer layers have been employed to ensure good con¬ 
tact between the different films. The island is finally 
Al(3nm)/Ti(2nm)/Au(10nm)/Ti(20nm). The resulting 
multilayer island is expected to play a role in enhanc¬ 
ing the interference in the normal side of the junction, 
but a systematic study would be required to make defi¬ 
nite statements. In the following we will treat the island 
as a single composite metal with a rectangular shape of 
length L = 2.73 /im, width W = 0.4 /im and thickness 
77 = 35 nm. The design maximizes the overlap area A 
of the junction, leading to A = 0.5 /im^ for a volume 
of V = 0.038 yum^. Figure [T]c depicts the measured sam¬ 
ple. This large area combined to the light oxidation (0.46 
mbar of O 2 during I min 15 s) results in two series junc¬ 
tions, which we expect to have nearly equal normal state 
resistance = 755 Under this assumption the ap¬ 
plied voltage is split equally between the junctions and 
^SINIs(^) = ^nis(^/2) . 

The current / through one junction is expressed as 
m = 4 y * g{e) [/(e - eV) - f(e + eC)], (1) 


where we assume quasi-equilibrium on both sides of the 
junction, neglect charge imbalance on the superconduct¬ 
ing side, and choose the distribution function in the 
normal island to be the equilibrium Fermi-Dirac one 
f{E) = I/(e^/^®^ + I). The energy dependent conduc¬ 
tance ^(e) has the form (see, e.g. Ref. [9]): 


g(e)=g^^^{e) + 


^(A-|e|)A^ E^{2e) 

A^ — 

[2IU(6)] 


A^ — 2e‘^usR% 


( 2 ) 


where us are the densities of states in the normal and 
superconducting leads respectively, W{e) = i\/A^ — 
for |e| < A, W{e) = sgn(e) \/e‘^ - ISA for |e| > A and 
^BTK(^) conductance of a junction connecting two 
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FIG. 1. (a) I-V characteristics measured in 4-probe config¬ 
uration. The sample is thermally anchored to the mixing 
chamber of a dilution refrigerator. The bath temperature 
Tbath = 60 mK is measured using a calibrated ruthenium ox¬ 
ide thermometer. The inset (b) shows a close up near zero 
bias, where the step due to the Andreev current is clearly vis¬ 
ible. The dashed line is a fit using = 37.5 pA and Te ~ 80 
mK. The top inset (c) is a Scanning Electron Micrograph 
(SEM) of the sample. 


bulk leads derived by Blonder, Tinkham and Klapwiik 

(BTK^ 


.btk 




2T^0{A - |e|)A2 
[T2e2 + (2-T„)2(A2-e2) 

2r„^(|e|-A)|e| 

T„|e| + (2-r„)Ve2-A2 


(3) 


Tn are the set of transmission probabilities of the con¬ 
ducting channels. They relate to the junction normal 
state resistance Rn hy the Landauer formula 1/Rn = 

The disorder enhanced Andreev reflections are produc¬ 
ing two corrections to the conductance oc S 7 v,s'(co’) ap¬ 
pearing in Eq. ([^. The functions S 7 v,s'(co’) are expressed 

as double integrals of the Cooperons over the 

junction area 


Siv,s'(w) ^ ^ J j 


(4) 


and the Cooperons themselves are the solutions of the 
diffusion equation 


{-iuj + - r>iv,s'V^) Cysi'^) = S(r - r'). (5) 


Here are dephasing times in normal and supercon¬ 
ducting leads respectively and -D]v,s are the correspond¬ 
ing diffusion constants. 

Taking a typical value of dephasing time, ^ 1 ns, 
and estimating the diffusion constant, D vpl/S ~ 2 


cm^/s (here 0.32 x 10^ m/s is Fermi velocity in 

titanium, and I ^ 3 nm as in Ref. HU, we estimate the 
effective Thouless energy of the device to be 

^Th ^ ^ ^ ^ 


This value is small as compared to typical bias voltages 
and temperatures, which allows us to simplify the ex¬ 
pression for the current 0 . Namely, since the func¬ 
tion SAr(2e) quickly decays for energies |e| > E’xh, while 
S 5 ( 2 kF(e)) varies slowly at energies |e| < A, we can make 
the approximations (see HU 


SAr(2e) 

^s{2W{e)) 


Gtt 

Y 




4 1 

AHs \J A2 — e2 ’ 


(7) 


where Rs is the thickness of the superconducting film. 
Substitutin g this result in Eq. 0 for eV < A, we arrive 
at the resulu^^^, 


/(l/)=/,p(b) + /ivV)+/5V), (8) 

where 

lYY = ^ / [/(e - eb) - /(e + eb)], 

pV pV 

In{V) =/i^tanh^^, Is{V) = Is / (9) 

The current amplitudes g read 

/A ^ 3^ h 2h 

^ 2 e^v^VRl^ ^ e^usAHsR%‘ ^ ’ 

Assuming ~ 7. 10^^ m“^ as the titanium den¬ 

sity of states, one gets Ig = 19 pA and = 7.25 pA 
with the parameters of our sample device. The experi¬ 
mental value of = 37.5 pA is obtained by fitting the 
iow bias region of the I-V characteristics at base tem¬ 
perature (dashed biack iine of Fig. [^. The discrepancy 
between the experimentai vaiue and its theoreticai pre¬ 
diction might be a consequence of the iayered isiand as 
the two vaiues wouid coincide if one considers an effective 
thickness of the normai isiand to be = 7.4 nm instead 
of 35 nm. 

At iow temperature {k^T A) and for weakiy trans¬ 
parent junctions (T^ <C 1), one can approximateiy ex¬ 
press zero bias conductance in the form 



is an effective Dynes parameter accounting for the sub¬ 
gap ieakage current. Environment assisted tunneiin^ 
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FIG. 2 . Zero bias resistance as a function of temperature. 
The fitting model, indicated by black lines, uses Eqs. (11) and 
(15). The parameters are Rn = 755 ^2, A = 221.7 /xel/, 7 = 
0.8 X 10“^, = 37.5 pA. One more parameter is the external 

heat current Qext accounting for the deviation observed at low 
temperature (see inset). 


may also contribute to the phenomenological parameter 
7 . Nevertheless, with the value of Ig = 19 pA estimated 
above, one gets elgRjsf/A « 0.65 x 10“^ which agrees 
well to 7 = 0.8 X 10“^ in the experiment. 

Figure depicts the zero bias resistance, extracted 
from the I- V characteristics by numerical differentiation. 
The aluminum gap A = 221.7/ieV is used as a fitting pa¬ 
rameter at temperatures near the aluminum transition, 
where the total conductance is dominated by the quasi¬ 
particle tunneling. The responsivity in this regime where 
ksT A but where quasiparticle conductance still dom¬ 
inates reads 


^^0 _ I A a//cbT 

dT TV 


(13) 


The NIS thermometer responsivity will eventually vanish 
due to the sub-gap leakage, as the NIS resistance reaches 
the limit set by the 7 parameter. 

The responsivity changes sign for the temperature 
k^To ^ A/ In ^ • For temperatures below Tq An¬ 

dreev current is dominating the zero bias conductance, 
leading to the positive responsivity 


dRo 2k^ 


(14) 


which remains constant for k^T > Tth- This phe¬ 
nomenon is also called reentrance effect, it has been pre¬ 
dicted in the context of NI S junc tion theoreticall^ F^F^ i 
and measured experimentall}^^^. 

In order to reproduce the data presented in Fig. [^theo¬ 
retically, we need to take into account that the electronic 
temperature saturates at around 85 mK when cooling the 
bath (Tbath) below this temperature. This can be incor¬ 
porated in the modeling by considering the normal island 
of the device to be a free electron gas in quasi-equilibrium 


so that its temperature is the solution of the following 
Heat Balance Equation (HBE): 


Qa + 2 Qqp(C Te) + = SV(T“ - Tb7h) • (15) 


Although Andreev current does not transport heat across 
the barrier, it has been shown that i t pro duces power 
Qa = In X V into the normal islancP^^. ( 3 qp(H,Tg) 
represents the heating (or cooling) due to quasi-particle 
transport through one of the junctions. One can neglect 
this term near zero bias, but it dominates the left-hand 
side of Eq. (15) as soon as eV approaches A. The term 


on the right-hand side is the electron-phonon coupling, 
with a = b and E ^ 1 nW/K^//im^ for most metalS 
The phonon temperature is assumed to follow the bath 
temperature Tbath so that assuming an external heat load 
Qext = 170 aW, the electronic temperature saturates 


around Tsat = y Qext/^V = 85 mK. 

Although not less than 5 parameters are needed to 
describe the zero bias resistance as a function of tem¬ 
perature fully, only a single one (J^) accounts for the 
constant responsivity in the range from Txh/^B to Tq. 
Furthermore, this parameter can be known beforehand, 
by measuring a single I-V characteristic. 

Let us now extract the responsivity and sensitivity of 
our implementation of an Andreev thermometer. Figure 
shows the voltage response of a current biased SINIS 
device. Parameters which have been extracted previously 
are re-used to plot the theoretical response (dashed lines) 
using Eqs. 0 and ( [T^ , excluding self heating. The up¬ 


per curve (/bias = 43 pA > /^) is a typical response of 
a NIS quasiparticle thermometer. The saturation, due 
to the 7 parameter, is limiting the range of this mea¬ 
surement, as suggested by Eq. 0 . On the contrary we 
do not expect saturation of the Andreev thermometer, 
whose maximum responsivity is given by 9K/9Tmax = 
2k^le 172 /iV/K. Due to the positive responsivity, 
implementing electro-thermal feedback requires a voltage 
biased device, but the relatively small logarithmic deriva¬ 
tive of the resistance, a = dinR/dhiT 1 will limit its 
strength. In order to reduce self heating one can operate 
the thermometer close to zero bias voltage and measure 
the impedance of the structure. The responsivity is then 
Iq = 4.28Mf2/K, which is in good agreement with the 
theoretical prediction given by Eq. (14). 


In the small bias regime, the main source of noise 
comes from the voltage pre-amplifier. The Noise Equiv¬ 
alent Temperature (NET) of the thermometer can then 
be directly calculated as 


NET = 


dV/dT 


= 20.8//K/\/Hz 


(16) 


where Sy is the input noise of the voltage amplifier, typ- 
ically « InV/v^. If one can assume the noise 
of the amplifier to depend only weakly on the sample 
impedance, the NET is expected to be constant over the 
full temperature range. 
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FIG. 3. Voltage response of the thermometer. The dashed 
lines are calculated using the parameters listed in Fig. [^cap¬ 
tion. The non-monotonic behavior of the voltage denotes two 
different regimes, quasiparticle NIS tunneling at high temper¬ 
ature and disorder enhanced Andreev current at low temper¬ 
ature. 


The Noise Equivalent Power (NEP) is related to the 
NET as NEP = Gth,E x NET^ where Gth,E is the total 
thermal conductance, linearized near the working point. 
Near zero bias voltage the electron-phonon coupling dom¬ 
inates the heat transport and Gth,E ~ bEVT^. We 
then expect NEP = 1 x 10“^^ at 100 mK, with 

E = 2.4 nW/K^//im^, which is the value of bulk gold. 
We believe this estimation of the NEP to be the worst 
case estimate as the island is composed partly of tita¬ 
nium with a smaller E value. As the responsivity of the 
Andreev thermometer is constant, reducing the temper¬ 
ature leads to a quick improvement of the NEP. At 20 
mK we would expect NEP = 1.6 x 10“^^ W/\/Hz within 
this model. 


Defining a figure of merit for a thermometer is a prob¬ 
lem with an application-specific answeil^. According to 
Eq. (14) the main advantage of Andreev thermometer 
appears to be its constant responsivity over a wide range 
of temperature, bounded below by the Thouless energy 
(Eq. §). A realistic estimation gives 2.5 ^ 25 mK de¬ 
pending on the value of = 1.5 ^ 0.15 ns measured 
in other experiment d^^ l ^^ i Once the Andreev current 
has been measured at a single bath temperature, the re¬ 
sponsivity is simply a number involving only fundamental 
constants. We measured the responsivity in our device to 
be dR/dT\Q = 4.28 MD/K, which was constant over the 
full temperature range of the measurement [80-200 mK]. 
Generally, in an equilibrium environment where Qext can 
be neglected, the saturation of the electronic temperature 
is determined by self heating. In this respect Andreev 


thermometer is favorable since it operates near zero bias. 
This is not the case for instance for a NIS thermometer 
for which a trade-off between the responsivity and the 
operation range cannot be avoided. 

As Andreev thermometer is probing the local temper¬ 
ature of a metallic island, one can expect this system to 
be a radiation absorber, and we estimate the NEP to be 
1 X 10“^^ W/V^Hz at lOOmK. Eurther measurements at 
lower temperature, by reducing the external heat load of 
the present experiment would be required to test the low 
temperature limitations of an Andreev thermometer. 

This material is based upon work supported by the 
Academy of Einland under projects no. 139172 and 
250280 (LTQ Centre of Excellence), and by the European 
Commission under project no. 264034 (Q-NET Marie 
Curie Initial Training Network). The research made use 
of OtaNano, the Otaniemi Research Infrastructure for 
Micro- and Nanotechnology. 

^Planck Collaboration I, Astron. Astrophys.571, A1 (2014). 

^J. Pekola, J. of Low Temp. Phys. 135, 723 (2004). 

^F. Giazotto, T. T. Heikkila, A. Luukanen, A. M. Savin, and J. 
P. Pekola, Rev. Mod. Phys. 78 , 217 (2006). 

^D. R. Schmidt, C. S. Yung and A. N. Cleland, Appl. Phys. Lett. 
83 , 1002 (2003). 

^S. Gasparinetti, K.L. Viisanen, O.-P. Saira, T. Faivre, M. Arzeo, 

M. Meschke, and J.P. Pekola, Phys. Rev. Appl. 3, 014007 (2015). 
®J. T. Muhonen, M. Meschke, and J.P. Pekola, Rep. Prog. Phys. 

75 , 046501 (2012). 

^J. P. Pekola, V. F. Maisi, S. Kafanov, N. Chekurov, A. Kemp- 
pinen, Yu. A. Pashkin, O.-P. Saira, M. Mottonen, and J. S. Tsai, 
Phys. Rev. Lett 105 , 026803 (2010). 

®S. Rajauria, P. Gandit, T. Fournier, F. W. J. Hekking, B. Pan- 
netier, and H. Courtois, Phys. Rev. Lett. 100 , 207002 (2008). 
^D. S. Golubev, M. S. Kalenkov, and A. D. Zaikin, Phys. Rev. 
Lett. 103 , 067006 (2009). 

^°G. E. Blonder, M. Tinkham, and T. M. Klapwijk , Phys. Rev. B 
25 , 4515 (1982). 

^^M. E. Gershenson, D. Gong, T. Sato, B. S. Karasik, and A. V. 

Sergeev, Appl. Phys. Lett. 79 , 2049 (2001). 

^^See supplemental material at [URL will be inserted by AIP] for 
the derivation of the Cooperon functions. 
i^F. W. J. Hekking and Yu. V. Nazarov, Phys. Rev. Lett. 71, 1625 
(1993). 

^"^F. W. J. Hekking and Yu. V. Nazarov, Phys. Rev. B 49, 6847 
(1994). 

15y, Tanaka, A. A. Golubov, and S. Kashiwaya, Phys. Rev. B 68, 
054513 (2003). 

^®Y.V. Nazarov and T. H. Stoof, Phys. Rev. Lett. 76 , 823 (1996). 
Quirion, C. Hoffmann, F. Lefloch, and M. Sanquer, Phys. Rev. 
B 65 , 100508 (2002). 

^®B. Reulet, A. A. Kozhevnikov, D. E. Prober, W. Belzig, and 
Yu.V. Nazarov, Phys. Rev. Lett. 90 , 066601 (2003). 

^®A. Bardas and D. Averin, Phys. Rev. B 52, 12873 (1995). 

^^A.S. Vasenko, E.V. Bezuglyi, H. Courtois, F.W.J. Hekking, Phys. 
Rev. B 81 , 094513 (2010). 

^^H. Pothier, S. Gueron, D. Esteve, and M. H. Devoret, Phys. Rev. 
Lett. 73 , 2488 (1994). 

^^F. Pierre, A. B. Gougam, A. Anthore, H. Pothier, D. Esteve, and 

N. O. Birge, Phys. Rev. B. 68, 085413 (2003). 












